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Abstract 

In crystalline color superconductivity, quark pairs form at non-zero total 
momentum. This crystalline order potentially enlarges the domain of color 
superconductivity in cold dense quark matter. We present a perturbative 
calculation of the parameters governing the crystalline phase and show that 
this is indeed the case. Nevertheless, the enhancement is modest, and to 
lowest order is independent of the strength of the color interaction. 
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While Quantum Chromodynamics (QCD) has now been firmly established as the theory 
of strong interactions, many important questions, such as the nature of confinement and 
chiral symmetry breaking, still remain. Attempts to address such issues have led to recent 
interest in exploring the nature of nuclear matter under extreme conditions of temperature 
and density. New phases of nuclear matter have long been sought after, both theoretically 
and experimentally. Under the right conditions, nuclear matter is expected to undergo a 
phase transition to a quark-gluon plasma. At cold temperatures, however, recent application 
of QCD suggests a much richer phase structure of quark matter than originally anticipated. 
Perhaps the most prominent of such phases, and one which may play an important role in 
the core of neutron stars, is that of color superconductivity 

Color superconductivity is essentially the quark analog of BCS superconductivity. Be- 
cause of attractive quark-quark interactions in QCD, the Fermi sphere of quarks becomes 
unstable against the formation of Cooper pairs and the system becomes superconducting 
at sufficiently low temperatures. Quantitative results may be obtained at ultra-high baryon 
densities where asymptotic freedom ensures the validity of the weak coupling approximation. 
In this regime, the main contribution arises from one-gluon exchange in the color antisym- 
metric channel. The resulting transition temperature, T c , is given in terms of the running 
QCD coupling g(ji) by J^O 



,9.5. / 6iVc -k~ z tt 2 +4, (at -|N 

k B T c = 5127TV 4 2 4 e ~ V ~~ C 1 ) 

where /x is the chemical potential and 7 = 0.5772 ... is the Euler constant. 

This result that T c ~ e- K l g is remarkable for its non-BCS behavior (where BCS theory 
instead predicts T^ BCS ^ ~ e~ x ^ 9 ). The relation between T c and the zero temperature gap 
energy Ao, on the other hand, remains BCS-like: Ao = ne'^ksTc. This non-BCS scaling 
is driven by the poor screening of the long range color magnetic interaction. The same 
mechanism is also responsible for non-Fermi liquid behavior which suppresses T c significantly 



Since the main attractive channel is between different species of quarks, search for a 
realistic color superconductivity phase must take into consideration the masses of the light 
quarks. This effect, not accounted for in Eq. ([!]), is especially important when considering 
the color flavor-locked state || involving u, d and s quarks. Although not strictly identical, 
this situation is well modeled by separating the Fermi surfaces of massless quarks that 
participate in the pairing. This case was recently investigated in Refs. ji2Hl5|, where it was 
shown that a new crystalline superconductivity state may arise for an appropriate window 
of separation between Fermi surfaces. 

Crystalline superconductivity was originally investigated in BCS theory, where ferromag- 



netic impurities separate the Fermi surfaces of each electron spin [jl6HT8|| , so that /ii 2 = fi±5. 
When the system is cooled from the normal phase, the transition temperature is the max- 
imum pairing temperature, allowing for a possible net momentum 2q for a pair. For a 
sufficiently small separation, 5 < 8 C , zero momentum pairing wins. However, for 5 in the 
range 5 C < 5 < 5 max , pairing with non-zero momentum is favored. The net momentum 
provides a region of phase space where both quarks may sit on their respective Fermi sur- 
faces. However it breaks rotational invariance and leads to crystalline order. This state is 
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also known as the LOFF state, after the work of Refs. []16ljL7| . Finally, for 5 > <5 max , the 
transition temperature drops to zero, and superconductivity is lost. 

In the BCS case, both 5 C and <5 max as well as the pair momentum g max (at 5 = <5 max ) 
are proportional to the initial (i.e. 5 = 0) zero temperature gap energy, with the con- 
stants of proportionality being independent of the details of the interaction. In particular, 
£(bcs) w 0.606 A , while <^ ax s) « 0.754 A . As this indicates, the window for BCS crystalline 
superconductivity is fairly narrow. 

In this letter, we explore the crystalline superconductivity state for color superconduc- 
tivity. Our approach follows the perturbative method developed in Refs. ||H which allows a 
systematic approach to determining T c (but not directly the zero-temperature gap). We find 
that for color superconductivity, while 5 C remains unchanged, <5 max is modestly enhanced, so 
that there is a widening of the region of crystalline order compared to the BCS case. As in 
the behavior of T c itself in Eq. (JJ), the resulting non-BCS behavior of the color LOFF state 
also traces its origin to the long range nature of the magnetic gluons. 

We begin with a brief review of the perturbative method developed in Refs. |||] for a 
normal-phase determination of the transition temperature, however now applied to pairing 
with total momentum 2q and shifted chemical potentials fi ± 5. The proper vertex function 
for quark pair scattering in the color antisymmetric channel is given by a Dyson-Schwinger 
equation 

T q ^(P f \P t ) = Tf/iPf^) + ^- 3 K,AP f \P)T,Am), (2) 

where the initial and final momenta are q±Pi and q±pf. Here P stands for Euclidean 
four- momentum, (p, — v) with v being the Matsubara energy, and the four-momentum of 
the pair is Q = (2q, 0). Eq. (H) is of the Fredholm type with kernel 

K^ S {P'\P) = T}y(P'\P)S(Q + P, ^ + 8)S(Q -P,ii-8), (3) 

where T 2PI denotes the two quark irreducible part of the vertex and S(P, /i) the full quark 
propagator with chemical potential \x. 

The pairing instability may be probed by examining the spectrum of the kernel, given 
by the eigenvalue equation 

Em = e/ 0^ K «AP\nf(n, (4) 

for eigenf unctions f{P)- At weak coupling, all eigenvalues are less than one for sufficiently 
high temperature. Then, as the temperature is lowered, a point is reached where the largest 
eigenvalue reaches unity, E max (T, fi,5,q) = 1, and the solution to Eq. @ becomes singular, 
corresponding to the onset of the pairing instability. 

We have indicated explicitly the dependence of -E max on the temperature, mean chemical 
potential, Fermi surface separation and the pair momentum. The condition -E max = 1 may 
thus be solved to obtain a pairing temperature as a function of q for given /i and 5; this 
provides a determination of T c (fi, 5) when maximized over q. Ordinary and crystalline order 
corresponds to q = and q ^ 0, respectively, at this maximum point. At weak coupling, the 
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relevant 8 and q are both much smaller than the mean chemical potential /x, and the onset 
of LOFF order can be treated as a perturbation from ordinary pairing at 8 = q = 0. Thus 
we look for an expansion of the eigenvalue E max in ascending powers of g (with coefficients 
that may depend on In g). Through Eq. ([!]), this is equivalent to an expansion in descending 
powers of In fi/ (/cbT c °), where T c ° is the 8 = value. Quantities such as k B T^ / fi, 8/fi and 
q/ fi are much smaller than any finite power of g and can be neglected. 

To the accuracy of Eq. (|I|), the contribution of antiparticles to the scattering amplitude 
can be ignored and the two quark irreducible part of the vertex Y 2PI is well approximated 
by single gluon exchange: 

r ?>s (P'\P) ~ -4(1+ ^)[D M (k,u)N M (p,p',k) + D E (k,u)N E (p,p',k)], (5) 

where k = p' — p and uo = v' — v. Here Nm and N E are kinematic factors and DM(k,oj) 
and D E (k,u) are magnetic and electric gluon propagators with hard dense loop (HDL) 
resummation. While the one- loop quark self energy contributes to the prefactor of Eq. (0), 
it will not interfere with the perturbative contribution of nonzero 8 and q. Consequently we 
replace S with the free quark propagator, S(u,p; fi) = %/{iv — p + fi), throughout. 

The eigenvalue problem, Eq. (|j), was previously analyzed within each angular momentum 
channel for 8 = q = 0. The highest eigenvalue with angular momentum J reads 

^ = 1 + A(3c. + ln^)+0(ln- 2 i), (6) 

e c 

where 

f i Pj(x) - 1 JO, if J = 0; 

and e = k B T°/jl where fi = 25Q7r 3 g- 5 (2/N f ) ' /x is a rescaled chemical potential. The 
corresponding eigenfunction is 

f IP)- 2 " ,in ( w]nl A Y '^ m 

where v = v/jl and Yjm{P) are spherical harmonics. Since the level spacing for different J is 
of the same order as the perturbation to be introduced, a degenerate perturbation method 
has to be employed. 

Nonzero 8 and q modify the pole structure of the quark pair propagator of Eq. (|3|) in a 
similar manner as in the BCS case. In addition, a nonzero 8 will shift the Debye mass in the 
HDL resummed denominators of the vertex, and the numerators of the one-gluon vertex, 
Nm and iV^, depends on the total momentum q. Both corrections to the vertex function 
are of the order of (8/fi) 2 or (q/fi) 2 and may be ignored. 

To diagonalize the kernel, Eq. (^]), with nonzero 8 and q to 0(ln _1 -), we evaluate K 
between the 0(1) eigenfunction f JM (P) and its adjoint 

7jm(P) = S(P)S(-P)f JM (P) ~ u2 f ^\ r (9) 
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and define 



Ej IM >,JM = "|pE fj'M'iP'^siPWjMiP). (10) 

Using the eigenvalue equation, Eq. (|J), for 5 = q = and to leading order in g, we find 

Ej, M >,JM = Sjj'Smm* + ^j(J'M'\h op \JM). (11) 

The operator h op is given by 3c op + v op , with c op diagonal in the angular momentum repre- 
sentation, (J'M'\c op \JM) = ScjSjj'Smm 1 , and v op diagonal in the coordinate (angle) repre- 



sentation, 



T ,A\ _ ,/l i 



= in ^ + - Re ^(2 ■ 2^ - ^ cos ^))' ( i2 ) 



where V(C) = ^ m ^(C)- Here 6 is the polar angle with respect to momentum q. The 
azimuthal quantum number M remains conserved. Both c op and Hermit ian and c m 

is non-positive. For q = 0, the rotation symmetry is intact and the largest eigenvalue of h op 
still corresponds to J = 0. Here we recover the BCS result 



E = 1 + 



2 



Inf + ^)-Re^-^)]. (13) 



With <f 7^ 0, different angular momentum components start to mix and this widens the 
LOFF window with respect to 5. In what follows, we shall focus on the component M = 
since the additional sign change introduced by nonzero M will not be the most favorable. On 
returning to a coordinate representation, Yjm(p) — > u(x = cos 9), the eigenvalue equation 



for h op becomes 



i u(x') — U{X) 



dx ; ; \~ V op (x)u(x) = Xu(x), (14) 

-1 \x — x'\ 

where the first term on the left hand side is the coordinate representation of c op , following 
from the integral representation of Cj in Eq. (|7|). 

To explore the upper threshold of the LOFF window, 5 max , we take the limit T — > 0, 
whereupon the function v op (x) reduces to 

M^)=^(^)-ln^ J -ln|l-^|, (15) 

with x c = 5/q. The only possibility for obtaining a diverging <5 max is to take full advantage 
of the logarithmic singularity of Consider, e.g., a trial function 

i \ J -kri for x c — a < x < x c + a; 
u(x) = < c (16) 

I 0, otherwise. 
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FIG. 1. The eigenfunction u(x) corresponding to the largest eigenvalue of Eq. (14). For fixed 
5, this occurs at x r = 0.863. 



In the limit a — > 0, we find v)v op u — > In - but v)c op u — > — In -, so that v)h op u — > —2 In -. 
Therefore the operator /i op is bounded from above, and the widening of the LOFF window 
cannot grow indefinitely. A numerical analysis reveals that the largest eigenvalue of the 
integral equation, Eq. flU]) , with the limiting v op of Eq. flloT) for a fixed <5 is A ~ — In 27rA .^ T o — 
1.304, which occurs at q ~ 1.16 5. The corresponding u(x) is shown in Fig. [I], demonstrating 
the peak at x — x c . For a second order phase transition, the eigenvalue equation coincides 
with the limit of the gap equation as T — > T c . An anisotropy will develop for the order 
parameter in the superphase as well. Applying the pairing condition, E max = 1, and the 
relation between T® and Ao, we find 

<5 max ~ 1.67fc fl T° ~ 0.96 A , (17) 

to be compared with the BCS case, ~ 0.754 A . 

The determination of the onset, 5 C , is numerically more tedious, but an upper bound 
can be readily obtained. For very small 5, we have v op ~ ln^o — -^r/3 2 (<5 — qcosO) 2 . The 
last term is too small to mix different J, and we have J = for pairing. Since (Qpr)s,p < 
at q = 0, the LOFF state is not favored. Then if we assume the onset of the LOFF state 
starts from q c = for sufficiently large 5, the critical value of S follows from the condition 
(f^)<5,/3 = 0. Applying Eq. (|13"D and the pairing condition E max = 1, we obtain 

6 C ~ 0.606 Aq, (18) 

which is identical to the BCS result. Although we have not been able to verify it analytically, 
numerical results indicate that pairing indeed starts at q = 0. The numerically computed 
phase diagram is shown in Fig. |2|. 
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We note that Ref. \\A\ suggested that there would be an indefinite widening of the 



window for the color LOFF state by drawing an analog with a 1 + 1 dimensional system on 
account of the forward singularity of the scattering vertex. The present results, however, 
demonstrate that while the window is indeed widened over that of the BCS case, the widening 
is fairly modest, as indicated in Fig. [2| This result may be understood since the reduction 
of the available phase space for pairing with increasing 6 in higher dimensions dominates 
over any enhancement arising from a 1 + 1 dimensional analog. Furthermore, unlike the 
isotropic BCS case, different angular momenta have to be considered simultaneously towards 
a determination of <5 max , and the long range order parameter itself would not be spherical. 
To our knowledge, such a mixing effect is entirely new and the approximations we have used 
are well controlled. 

This angular momentum mixing in the color LOFF order is quite unique to the dynami- 
cally screened magnetic interaction. On writing the highest eigenvalue of the kernel, Eq. @, 
with angular momentum J, Ej = + 0(ln _1 -), we have Ef^ = 1 for all J here. On the 
other hand, for a BCS model with a pointlike interaction, pairing only occurs in the s-wave 
channel, and Ef^ = 5,;^. Even for a Yukawa type of interaction, say the screened electric 
interaction alone, E^ would be different for different J. In both the pointlike and Yukawa 
cases, the angular momentum only mixes in higher order corrections, and would not lead to 
the enhancement seen in the color LOFF state. 



7 



ACKNOWLEDGMENTS 

We wish to thank K. Rajagopal for providing useful insight and comments. This work 
is supported in part by the US Department of Energy under grants DE-FG02-91ER40651- 
TASKB and DE-FG02-95ER40899-TASKG. 



8 



REFERENCES 



[1] B. Barrois, Nucl. Phys. B 129 (1977) 390; S. Frautschi, Proceedings of the Workshop on 
Hadronic Matter at Extreme Energy Density, N. Cabibbo, ed., Erice, Italy (1978); D. 
Bailin and A. Love, Phys. Rept. 107 (1984) 325, and references therein for early works. 

[2] M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537 (1999) 443; M. Alford, 
|hc-p-ph/0 1020471 ; K. Rajagopal and Wilczek, [hep-ph/0011333| , to appear in B.L. Ioffe 
Festschrift, At the Frontier of Particle Physics/Handbook of QCD, M. Shifman ed., 
(World Scientific 2001), and references therein. 

[3] D.T. Son, Phys. Rev. D 59 (1999) 094019. 

[4] T. Schafer and F. Wilczek, Phys. Rev. D 60 (1999) 114033. 

[5] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61 (2000) 051501. 

[6] R.D. Pisarski and D.H. Rischke, Phys. Rev. D 61 (2000) 074017. 

[7] D.K. Hong, Nucl. Phys. B 582 (2000) 451. 

[8] W. Brown, J.T. Liu and H.C. Ren, Phys. Rev. D 61 (2000) 114012. 

[9] W. Brown, J.T. Liu and H.C. Ren, Phys. Rev. D 62 (2000) 054016. 
[10] W. Brown, J.T. Liu and H.C. Ren, Phys. Rev. D 62 (2000) 054013. 
[11] Q. Wang and D.H. Rischke, Phys. Rev. D 65 (2002) 054005. 
[12] M. Alford, J. Bowers and K. Rajagopal, Phys. Rev. D 63 (2001) 074016. 
[13] J. Bowers, J. Kundu, K. Rajagopal and E. Shuster, Phys. Rev. D 64 (2001) 014024. 
[14] A. Leibovich, K. Rajagopal and E. Shuster, Phys. Rev. D 64 (2001) 014024. 
[15] D. Hong and Y. Sohn, Critical Couplings in Crystalline Color Superconductivity, [hep 



ph/0107003| ; D. K. Hong, Crystalline Color Superconductivity in Dense Matter, |hcp- 
ph/0112028| . 



[16] A.I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz 47 (1964) 1136 [Soviet Physics 

JETP 20 (1965) 762]. 
[17] P. Fulde and R. A. Ferrell, Phys. Rev. 135 (1964) A550. 
[18] S. Takada and T. Izuyama, Prog. Theor. Phys. 41 (1969) 635. 



9 



